Mir^ncOPV  9^1 1 IT N  I FS1  CHART 

i‘.  I'th'  * 


I  -4 ■>**  j 


t'JU- ■Jj 


OTIC  FILE  CJti 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  w  *  »U  I  1 1.4.  SjU  Pi 


REPORT  DOCUMENTATION  PAGE 


VI  \T  I  lb  RESTRICTIVE  MARKINGS 


AD-A188  146 


)TIO 


mi 2 


3  DISTRIBUTION /AVAILABILITY  OR  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

USAETL 

6b.  OFFICE  SYMBOL 
(If  applicable) 

CEETL-LO) 

7a  NAME  OF  MONITORING  ORGANIZATION 

6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

7b  ADDRESS  (City,  State,  and  ZIP  Code) 

Fort  Belvoir,  VA  22060-5546 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

8b  OFFICE  SYMBOL 
(if  applicable) 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

1  8c.  ADDRESS  (City,  State,  and  ZIP  Code)  1 

10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 

PROJECT 

TASK 

ELEMENT  NO 

NO. 

NO 

WORK  UNIT 
ACCESSION  NO. 


1 1  title  (Include  Security  Classification)  Estimation  of  Gravity  Vector  Components  from  Bell  Gravity 

Gradiometer  and  Auxiliary  Data  under  Consideration  of  Topography  and  Associated  Analytical 
Upward  Continuation  Aspects. 


12.  PERSONAL  AUTHOR(S) 

Dr.  H.  Baussus  von  Luetzow 


13a.  TYPE  OP  REPORT  IT 3b  TIME  COVEREO  14  DATE  OF  REPORT  (Year,  Month,  Day)  IIS  PAGE  COUNT 

Research  I  FROM _ TO _  29  JAN  87  I 


17. _ COSATI  CODES _  18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

FIELD  I  GROUP  I  SUB-GROUP  .  _  ..  _  .. 


^  gravity  gradiometer  applications 
"topographic  noise" 


19.  ABSTRACT  ( Continue  on  reverse  if  necesfary  and  identify  by  block  number ) 

Following  an  introduciton,  the  paper  briefly  discusses  gravity  gradiometer  applications. 

It  outlines  the  estimation  of  first  order  derivatives  of  the  anomalous  gravity  potential 
from  Bell  gravity  gradiometer  and  auxiliary  data  in  the  context  of  a  Wiener-Ko lmogorov 
optimization  scheme  under  consideration  of  computable  "topographic  noise,"  accomplished  on 
the  basis  of  the  Pellinen-Moritz  solution  of  the  boundary  value  problem  of  physical 
geodesy.  The  paper  also  addresses  four  different  methods  of  analytical  upward  continuatior 
of  first  order  derivatives  of  the  anomalous  gravity  potential  under  identification  of  a 
finite  difference  method  using  Laplace's  equation  as  the  most  economical  and  efficient 
one.  Relevant  conclusions  are  then  presented. 


20.  DISTRIBUTION  /  AVAILABILITY  OF  ABSTRACT 

□  unclassified/unlimited  EFsameasrpt 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
F.  DARLENE  SEYLER 


□  OTIC  USERS 


ITY  CLASSIFICATION 
D 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 
202-355-2647  CEETL-LO 


CEETL-LO 


00  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Estimation  of  Gravity  Vector  Components  from  Bell  Gravity 
Gradiometer  and  Auxiliary  Data  under  Consideration  of  Topography 
and  Associated  Analytical  Upward  Continuation  Aspects. 


H.  Baussus  von  Luetzow 

U.S.  Army  Engineer  Topographic  Laboratories 
Fort  Belvoir,  VA  22060-5546 

Abstract 

Following  an  introduction,  the  paper  briefly  discusses  gravity  gradiometer 
applications.  It  outlines  the  estimation  of  first  order  derivatives  of  the 
anomalous  gravity  potential  from  Bell  gravity  gradiometer  and  auxiliary  data  in  the 
context  of  a  Wiener-Kolmogorov  optimization  scheme  under  consideration  of  computable 
"topographic  noise,"  accomplished  on  the  basis  of  the  Pellinen-Moritz  solution 
of  the  boundary  value  problem  of  physical  geodesy.  The  paper  also  addresses 
four  different  methods  of  analytical  upward  continuation  of  first  order  derivatives 
of  the  anomalous  gravity  potential  under  identification  of  a  finite  difference 
method  using  Laplace's  equation  as  the  most  economical  and  efficient  one. 

Relevant  conclusions  are  then  presented. 


1 ;  ,  INTRODUCTION 


The  estimation  of  gravity  vector  components  from  Bell  gravity  gradiometer 
and  auxiliary  data  has  been  discussed  by  a  number  of  authors,  including  Jordan 
(1982),  White  and  Goldstein  (1984),  Center,  Jordan  and  Peacock  (1985),  and 
Baussus  von  Luetzow  (1985).  Apart  from  gravity  gradiometer  self  noise  power 
spectra  or  corresponding  covariance  functions,  the  collocation  theory  of 
physical  geodesy  was  employed  to  derive  spectra  and  cross-spectra  or  equivalent 
covariance  functions  involving  first  and  second  order  derivatives  .of  the  anomalous 
gravity  potential.  The  basic  covariance  functions,  pertaining  either  to  a 
randomized  anomalous  gravity^  potential  or  to  a  randomized  gravity  anomaly, 
accomplished  by  the  subtraction  of  low’ frequency  harmonic  representations  of 
anomalous  potential  or  gravity  anomaly  data,  presuppose  the  existence  of  homogeneity 
and  isotropy.  The  application  of  collocation  theory  does  not  face  difficulties 
in  the  case  of  quasi-flat  terrain  and  is  a  prerequisite  to 
arrive  at  small  estimation  errors  for  gravity  vector  components 
in  the  context  of  an  airborne  survey  over  a  square  area  of  300  x  300  km* ,  covered 
with  parallel  traverses  and  cross  traverses  5  km  apart,  conducted  at  an  elevation 
of  600  m  at  a  speed  of  about  360  km  hr  ,  and  employing  only  a  limited  amount  of 
measured  gravity  vector  component  data.  It  is,  however,  well  known  that  the 
assumption  of  homogeneity  and  isotropy  does  not  satisfactorily  hold  in  the  presence 
of  moderate  to  strong  mountainous  terrain.  Some  efforts  have  been  made  to  apply 
heuristic  topographic  corrections  to  gravity  gradiometer  measurements.  Still, 
the  problem  remains  to  provide  for  a  Wiener-Kolmogorov-type  estimation  process 
free  of  computable  topographic  noise  so  that  the  assumption  of  homogeneity  and 
isotropy  is  reasonably  fulfilled  with  respect  to  modified,  signal-type  measurements 
disturbed  by  gradiometer  self  noise  only.  The  solution  to  this  problem  is 
possible  under  application  of  approaches  to  the  geodetic  boundary  value  problem 
and  the  correlated  interpolation  of  gravity  anomalies  and  deflections  of  the 
vertical  in  mountainous  terrain  by  Pellinen  (1969),  Moritz  (1969),  and  Baussus 
von  Luetzow  (1971,  1981).  Section  3  as  the  main  section  of  this  paper  addresses 
the  solution  under  consideration  of  the  estimation  process  and  the  computation 
of  "topographic  noise".  The  preceding  section  addresses  in  a  compact  way 
gravity  gradiometer  applications  to  place  the  determination  of  gravity  vector 
components  over  relatively  large  areas  and  the  analytical  upward  continuation 
of  vertical  and  horizontal  derivatives  of  the  anomalous  gravity  potential  in 
the  proper  perspective.  The  relatively  short  sections  4-7  discuss  different 
methods  of  analytical  upward  continuation,  followed  by  several  conclusions. 

2 .  GRAVITY  GRADIOMETER  APPLICATIONS 


Geodetic  research  and  the  development  of  associated  technology  has  been 
considerably  stimulated  by  military  requirements  and  actual  or  potential 
military  applications  with  critical  funding  provided  by  the  U.S.  Government. 

The  result  thereof  was  an  increase  of  scientific  knowledge,  new  instrumentation 
and  techniques,  and  technology  transfer  to  the  civilian  sector,  particularly  in 
navigation  and  surveying. 

The  main  applications  of  gravity  gradiameters  follow; 

•  Establishment  of  the  spatial  gravity  field  in  ICBM*  launch  areas  (most  important 
application).  This  includes  gravity  vector  components  on  the  ground  for 
calibration,  local  to  global  coordinate  transformations,  geodetic  azimuth 
determination,  and  gravity  programmed  inertial  positioning  (beneficial  to 
Small  ICBM  or  Midgetman) .  An  accurate  spatial  gravity  field,  extending  to 
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Intercontinental  Ballistic  Missile 


a  height  of  200  km,  provides  also  for  comparisons/calibration  frr  the 
Geopotential  Research  Mission  (GRM)  and  a  Superconducting  Gravity  Gradiometer 
Mission  (SGGM).  Instrument:  Bell  Gravity  Gradiometer. 
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0  Accuracy  enhancement  of  SLBM's  .  Instrument:  Bell  Gravity  Gradiometer. 

0  Improvement  of  global  gravity  field  for  geophysical  applications  (GRM). 
Instrument:  Superconducting  Gravity  Gradiometer. 

Byproducts:  ICBM  target  error  reduction  (not  critical),  cruise  missile 
navigation,  improved  satellite  ephemeris. 

0  High  accuracy  space  navigation.  Instrument:  Superconducting  Gravity 
Gradiometer  (six  axis  version). 

°  Gravity  vector  densif ication  over  large  areas  for  gravity  programmed 

inertial  navigation  and  for  geodetic  network  adjustments.  Instrument:  Bell 
Gravity  Gradiometer. 

°  Subterraneous  mass  detection  (cavaties,  oil,  water).  Instrument:  Super¬ 
conducting  Gravity  Gradiometer  (vertical  axis  version). 


°  Test  of  Newton's  square  law  and  tests  of  Einstein's  theory  of  relativity. 
Instrument:  Superconducting  Gravity  Gradiometer. 

3.  DETERMINATION  OF  DERIVATIVES  OF  THE  ANOMALOUS  GRAVITY  POTENTIAL  OVER 
MOUNTAINOUS  TERRAIN  USING  BELL  GRAVITY  GRADIOMETER  MEASUREMENTS  AND 
AUXILIARY  DATA 


A.  Notations  and  Relations  Used 


anomalous  gravity  potential 

derivatives  taken  along  the  local  horizon  in  a  northern  and 
eastern  direction 


fl  *  **,  Bell  gravity  gradiometer  instrument  noise  with  white  noise 
component  n >  and  red  noise  component 


A 


signal  variable  indicator 


% 

SI' 
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topographic  noise  component  of  Tx 
topographic  noise  component  of 

lover  order  spherical  harmonics  representation  of  7~x  used  for 
subtraction 


lower  order  spherical  harmonics  representation  of  7T  used  for 
subtraction 


lower  order  spherical  harmonics  representation  of  17*  used 
for  subtraction 


Sea  Launched  Ballistic  Missile 
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lower  order  spherical  harmonics  representation  of 
used  for  subtraction 


Z'SZ 

z-z  =  z  -*z-(z-st.) 
r„-sr„ 

Z»~  Z,  -  4*  ~^Zx ~(z*~  fZ* ) 


weight  factor  with  measurement  index  t  and  computation 
index  K  ,  applicable  to  7^  in  Wiener-Kolmogorov  estimation 

weight  factor  with  measurement  index  J  and  computation 
index  K  ,  applicable  to  f  in  Wiener-Kolmogorov  estimation 

earth's  mean  radius 


unit  sphere  (full  solid  angle) 


length  of  radius  vector  from  earth's  center  to  a  moving  point 
on  earth's  surface 

length  of  radius  vector  from  earth's  center  to  a  fixed  point 
on  earth's  surface 

length  of  radius  vector  from  earth's  center  to  a  fixed  point 
on  a  level  surface  above  earth's  surface 

— »  —T 

angle  between  radius  vectors  r  and  K  or  r 

f  A 

(  rz  y  rA  -  2r  r  cmP 

A  * 

azimuth  angle  counted  clockwise  from  north 


S(r \,*r)  generalized  Stokes  Function 


elevation  of  terrain  referring  to  a  movable  point  on  earth's 
surface 

elevation  of  terrain  referring  to  a  fixed  point  on  earth's  surface 
elevation  of  level  surface  above  terrain 


<  =  Z  %<:»!£ 

a  Z 


e,-  [frHJ]  'h 


planar  approximation  to  £ 


gravity  anomaly  on  earth's  surface 
topographic  correction  to  gravity  anomaly 
Faye  anomaly 


Bouguer  gradient 
gravitational  constant 
standard  density 


ST  -  + 


potential  of  the  topographic  masses 


ti-i. ) 


with  V  as  an  arbitrary  variable 


<P 

A 

3  -i-t 


geographic  latitude 


geographic  longitude 


B.  Estimation  Process 


A  tfiener-Kolmogorov  estimation  can,  under  consideration  of  /A  and  J" 
only,  be  characterized  by  ** 


Kit  "1  r*‘  *  %  iv  0«j  +  nj  J 


iVvyvVvppfi9v?i^p*ra*v^ivT^v*juii"iii  'um  — - 


The  derivatives  i k l  ,  also  called  truth  data,  are,  of  course, also  associated 
with  measurement  errors.  With  respect  to  deflections  of  the  vertical, 
proportional  to  /r  and  Ty  ,  efforts  are  being  made  to  determine  them 
astrogeodetically  with  an  error  of  0.1  arcsec  rms.  Gravity  anomalies 
can  be  measured  with  a  relatively  greater  accuracy.  In  order 

to  compute  the  vertical  derivative  ££  *  Tr  ,  high  degree  spherical  harmonic 
expansions  presented  by  Rapp  (1987)  or  the  GR1M3-L1  Model  described  by  Reigber 
Balmino,  Muller,  Basch  and  Moynot  (1985)  may  be  employed.  These  models  are 
also  useful. for  the  subtraction  of  lower  harmonics  components  from  measured 
first  and  second  order  derivatives  of  7"  . 

Regarding  eq.  (1),  it  is  assumed  that  all  measured  quantities  refer  to  a 
1*V®1  surface.  In  the  case  of  quasi~flat  terrain,  “estimates  have  been 

made  from  on  the  earth's  surface  and  from  airborne  Tr*  “measurements, 
employing  spatial  covariance  functions.  This  is  not  possible  in  the  presence 
of  pronounced  topography.  In  this  case,  it  is 


a) 


where  -pffi 
measurements . 
used. 


is  a  representative  gradient  obtained  from  gravity  gradiometer 
There  is  no  significant  degradation  in  accuracy  if  eq.  (2)  is 


For  simplicity,  eq.  (1)  presupposes  the  prior  application  of  quasi-systematic 
corrections  obtained  from  lower  harmonic  expansions.  In  the  presence  of 
pronounced  mountainous  terrain,  eq.  (1)  is  reformulated  as 


C.  Computation  of  "Topographic  Noise" 

The  essential  task  is  the  computation  of  topographic  noise  components 
of  the  derivatives  7r  ,  7y  ,  7%  •  In. this  respect,  the  solution  nomenclature 
employed  by  Moritz  (1969)  is  applied.  In  order  to  obtain  solutions  at 
elevation  *>A  instead  of  solutions  at  elevation  6  *  hp  ,  Stokes  generalized 
function  Oj  replaces  Stokes  function  SC'f)  in  the  solutions  where  applicable.4 

First,  the  solution  for  T  at  a  point  P^  is  written  as 


/ 


3 

4 


Moritz  (1969),  p.  31. 

Sj  i»  explicitly  formulated  by  Moritz  (1966). 


See  p.  49. 


The  solution  components  are 


Tf  represents  the  "signal"  term  of  /  since  its  structure  is  compatible 
with  the  existence  of  a  reproducible  kernel  in  the  sense  of  Krarup,  a 
prerequisite  for  the  application  of  the  collocation  method  of  physical  geodesy. 

and  are  clearly  "noise"  terms  which  do  not  permit  the  establishment 
of  a  statistical  estimation  structure.  ^  appears  to  be  a  more  slowly  varying 
function.  Moritz  (1969)  states  that  the  te_rm  ST  may  be  neglected  in  the  context 
of  a  planar  approximation.  Accordingly,  /j  appears  to  be  mainly  responsible 
for  fluctuations  of  a  non-stationary  nature. 

An  improvement  of  the  above  analysis  is  possible  by  the  replacement  of 
by  the  isostatic  anomaly  ,  provided  that  isostatic  conditions  prevail. 

Neglecting  a  second  order  improvement  involving  consideration  of  6f  -terms,  the 
additional  "topographic  noise"  can  be  written  as 

~  vr  Jj  ^a3f  ~  ^  to) 

* 

Without  consideration  of  the  term  ^  ,  the  "topographic  noise"  terms 
involving  vertical  and  horizontal  derivatives  of  7)  and  7^  can  be  formulated  as 
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Moritz  (1969),  p.  32. 
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D.  Averaged  Second  Order  Derivatives 


Because  of  the  application  of  a  moving  average  of  gravity  gradiometer 
measurements  over  a  time  interval  a/  of  about  10  seconds,  it  is  not  necessary 
to  compute  second  order  derivatives  with  respect  to*_.X  *  Z  of  7J"  and  7$  •  As 
an  example,  the  averaged  second  order  derivative  of  /j  with  respect  to  X  is 


In  eq.  (13)  it  is  A  X  m  ^  ~ J  » 

the  horizontal  aircraft  speed. 


&s) 


~tS^  being 


4.  STANDARD  ANALYTICAL  UPWARD  CONTINUATION 


Assuming  that  the  three  first  order  derivatives  of  the  anomalous  gravity 
potential  have  been  estimated  from  gradiometric  surveys  at  rectangular  grid 
points  of  a  level  surface  at  elevation  It ^  -const.,  Poisson's  integral  formula 


can  be  used  for  the  analytical  upward  continuation  of  T  or  of  its  derivatives 

*  4*  •  According  to  Heiskanen  and  Moritz  (1967),  for  elevations 

smaller ’than  250  km,  the  planar  form  of  eq.  (14)  is  sufficiently  accurate.  For 
elevations  up  to  150  km,  which  apply  in  ICBM  launch  areas,  the  horizontal 
integration  distance  has  to  be  ten  times  the  elevation  to  assure  good  accuracy, 
i.e.,  a  very  large  information  base  is  required.  Additionally,  rectrangular 
grid  data  has  to  be  interpolated  to  conform  with  circular  zone  data  to  be 
employed  in  the  computation. 

5.  IMMEDIATE  SPATIAL  WltNER  KOtMOGOROV  ESTIMATION 


The  three  derivatives  of  the  anomalous  gravity  potential  may  be  estimated 
directly  for  grid  points  at  higher  elevations  in  the  context  of  a  Wiener-Ko ltnogorov 
scheme,  using  spatial  covariance  functions  or  equivalent  power  spectra.  For 
higher  elevations,  data  derived  from  a  greater  number  of  adjacent  surveys  would 
be  required,  and  matrix  inversions  would  become  very  complicated.  Further 
complications  would  arise  in  the  presence  of  pronounced  topography  because 
of  the  non-applicability  of  effective  collocation  methods. 

6.  SPATIAL  COLLOCATION  UPWARD  CONTINUATION  OF  4.  .  /<*  .  71 


if  %  ,rx  are  first  estimated  at  grid  points  of  a  level  surface  with 
elevation  ,  a  collocation  upward  continuation  would  be  a  replacement  of 
the  Poisson  integral' formula  approach  by  statistical  methods.  Again,  a  large 
information  base  would  be  required,  inversions  of  matrices  associated  with 
higher  elevations  would  become  complicated,  and  complications  would  exist  in 
the  presence  of  pronounced  topography.  It  would  be  possible  though  to  apply 
the  collocation  extrapolation  only  to  signal  quantities  /r  ,  ,  A  and 

to  compute  the  "topographic  noise"  components  separately  for  higher  elevations 
by  means  of  eqs.  (9)  -  (12). 

7.  NUMERICAL  UPWARD  CONTINUATION  USING  LAPLACE’S  DIFFERENTIAL  EQUATION 


The  upward  continuation  methods  discussed  in  sections  4-6  are  characterized 
by  the  use  of  a  specific  derivative  of  the  anomalous  gravity  potential,  given 
at  a  level  surface,  by  the  requirement  of  a  large  information  base  for 
extrapolations  to  higher  elevations,  by  associated  large  inversion  matrices 
in  collocation  approaches,  and  by  the  inadequacy  of  these  in  the  presence  of 
mountainous  terrain  in  the  absence  of  a  separation  into  signal  estimation 
and  "topographic  noise"  computation.  Fortunately,  the  availability  of  three 
vector  components  at  a  level  surface  makes  it  possible  to  use  information  at 
two  adjacent  level  surfaces  as  a  prerequisite  for  applying  Laplace's  equation. 
Hereby,  the  information  base  can  be  reduced  significantly,  and  inversions 
of  matrices  are  not  needed. 


In' Cartesian  coordinates,  the  avai 
s  at  the  levelZg  makes  it  possible 
the  level  s  Ze  t  AZ  since 

jSL  kZ 

ZZ  dx  J*  Jz 


liability  tt  %  .  % ,  jT.  T 
sible  to  compute  T  .  /,  »  A 


2at  A  Z 


since 


M.  U  .  A  U 

ZZ  dX  d*  JZ 

&  % 
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UU) 

(/ft) 


The  derivative  "fijA  can»  of  course,  be  determined  separately  through  use  of 
gradiooeter  information,  and  a  representative  value  can  be  employed  for 
analytical  upward  continuation. 

If  a  variable  -f  satisfying  Laplace's  equation  is  given  ftjt  two  adjacent 
levels  separated  by  aX  ,  Laplace's  equation 


permits  application  of  the  numerical  algorithm 


6  -  *f.  -/•  -  n> 

In  more  appropriate  spherical  coordinates  T  ,  Q  ,  A.  ,  eq.  (17)  has  to 
be  replaced  by 

f*  =  lf'  '  /*  (2*¥r  t  c‘^6¥>  (/t) 


where  the  derivatives  in  the  bracketed  term  are  to  be  evaluated  by  finite 
differences.  Higher  order  finite  difference  algorithms  may  be  applied  if 
considered  advantageous. 

It  is  possible  to  use  eq.  (18)  only  for  "signal"  components  and  to 
compute  the  "topographic  noise"  effects  separately,  using  grid  lengths  of  1  km 
instead  of  5  km. 

8.  CONCLUSIONS 

Wiener-Kolmogorov-type  estimation  of  gravity  vector  components  or  of 
related  first  order  derivatives  of  the  anomalous  gravity  potential  from  Bell 
gravity  gradiomcter  and  auxiliary  data  in  the  context  of  an  airborne  area 
survey  can  also  be  accomplished  in  the  presence  of  pronounced  mountainous 
terrain.  A  specific  solution  of  the  boundary  value  problem  of  physical 
geodesy  permits  the  determination  of  "signal"  and  "topographic  noise" 
solution  terms.  "Topographic  noise"  terms  can  be  computed  for  each  first 
derivative  of  the  anomalous  gravity  potential,  and  their  consideration  permits 
the  application  of  algorithms,  including  numerical  weight  factors,  valid  in  the 
case  of  quasi-flat  terrain.  A  slight  degradation  in  accuracy  may  be  expected. 
The  standard  analytical  upward  continuation  method  and  two  spatial  collocation 
methods  for  the  estimation  of  first  order  derivatives  of  the  anomalous  gravity 
potential  can  be  replaced  by  numerical  upward  continuation  using  Laplace's 
equation  because  of  the  availability  of  three  gravity  vector  components  at 
the  information  base  level  surface.  In  this  case,  the  standard  information 
base  can  be  reduced  considerably,  and  no  matrix  inversions 

are  required.  Analytical  upward  continuation  of  "topographic  noise"  components 
may  also  be  accomplished  separately,  using  a  higher  resolution  grid.  The 
computation  of  "topographic  noise"  effects  is  mathematically  laborious,  can; 


however, be  achieved  by  means  of  high  speed  computers.  The  alternative  in  the 
presence  cf  mountainous  terrain  is  to  reduce  flight  traverses  from  300  km  to 
about  100  km,  to  use  highly  accurate  "truth"  data  at  both  ends  of  the 
traverse,  and  to  repeat  single  traverse  surveys. 
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